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We report on a search for B® — > /i + fjT and B® — > /i + fx~ decays in pp collisions at y/s = 1.96 TeV 
using 364 pb -1 of data collected by the CDF II detector at the Fermilab Tevatron Collider. After 
applying all selection requirements, we observe no candidates inside the B® or BjJ mass windows. 
The resulting upper limits on the branching fractions are B{B° -f < 1-5 x 1(T 7 and B(B° d -> 

M + /0 < 3.9 x 10~ 8 at 90% confidence level. 

PACS numbers: 13.20.He, 12.15.Mm, 12.60.Jv, 13.85. Qk 



In the standard model (SM), Flavor Changing Neu- 
tral Current (FCNC) decays are highly suppressed and 
can only occur through higher order diagrams. The de- 
cay rate for the FCNC decay B Q S — > 1} is propor- 
tional to the CKM matrix element | T/f s | 2 . The rate of 
B® — > decays is further suppressed by the ratio 
of |Vtd/Vt s | 2 . The SM expectations for these branching 
fractions are B{B° -> ^ + aO = (3-42 ± 0.54) x 1CT 9 and 
B(B° -► v+fi-) = (1.00 ±0.14) x 10~ 10 2], which are 
about two orders of magnitude smaller than the current 



experimental sensitivity. However, new physics contribu- 
tions can significantly enhance these branching fractions. 
An observation of these decays at the Tevatron would be 
unambiguous evidence for physics beyond the SM. The 
best existing experimental bound on B is < 4.1 x 10~ 7 Q 
(< 8.3 x 10~ 8 4]) for B° S (B°) -> fj+fi' at the 90% con- 
fidence level (C.L.). 

In minimal supersymmetric (SUSY) extensions of the 
SM, additional diagrams involving SUSY particles also 
contribute to FCNC decay rates and the branching frac- 
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tion B(B® d — > /i + /i~) oc (tan/3) 6 , where tan/3 is the 
ratio of vacuum expectation values of the two neutral 
CP-even Higgs fields. Large values of tan /3 enhance the 
decay rate to a level observable by the Tevatron experi- 
ments . For example, increases of one to three orders 
of magnitude are obtained in the minimal SO(10) mod- 
els , which favor large values of tan /3. For the minimal 
flavor violating (MFV) models, B d — > remains sup- 

pressed relative to B° — » due to |V t d/V ts | 2 . This 

may not be true for non-MFV models such as R-parity 
violating SUSY 7], which can produce large enhance- 
ments, even for low values of tan/3, in either or both of 
the B s and Bd FCNC decay rates. Thus, a simultaneous 
observation of B® d — > decays can be important in 

determining the flavor structure of the new physics. In 
the absence of an observation, any improvements to the 
limits can be used to constrain significantly many SUSY 
models 00011. 

In this Letter, we report on a search for B® — > 
and B d — » /i + /i - decays using 364 pb _1 of data collected 
by the upgraded Collider Detector at Fermilab (CDF II) . 
This data set includes 171 pb" 1 of data from our previous 
measurement 0. We significantly improve the sensitiv- 
ity of the search over our former analysis by doubling our 
data sample, extending the muon acceptance, and using a 
likelihood ratio technique for signal and background sep- 
aration. The limits we present here are the most stringent 
to date and supersede our previous results. 

The CDF II detector is described in detail in Ref. • 
The inner tracking system is composed of a silicon mi- 
crostrip detector (SVX II) 111 surrounded by an open- 
cell wire drift chamber (COT) [12j . These tracking detec- 
tors are immersed in a 1.4 T magnetic field and measure 
Pt, charged particle momentum in the plane transverse 
to the beamline. Four layers of planar drift chambers 
(CMU) detect muon candidates withpr > 1.4 GeV/c 
and provide coverage in the pseudorapidity range |ry| < 
0.6, where r\ = — ln(tan |) and 9 is the angle of the track 
with respect to the beamline. The central muon exten- 
sion (CMX) consists of conical sections of drift tubes and 
extends the coverage to 0.6 < \r/\ < 1.0 for muon candi- 
dates with p T > 2.0 GeV/c. 

The data used in this analysis arc selected by two 
classes of dimuon triggers: for the CMU-CMU (U-U) 
triggers both muon candidates are reconstructed in the 
CMU chambers, while for the CMU-CMX (U-X) triggers 
one of the muon candidates is reconstructed in the CMX 
chambers. The inclusion of the U-X trigger increases the 
signal acceptance by about 50%. The details of the trig- 
ger system and selection requirements can be found in 
Refs. 0E3- Since they have different sensitivities, we 
treat U-U and U-X channels separately, combining the 
results at the end. 

The offline reconstruction begins by identifying two 
muon candidates of opposite charge which satisfy the 
online dimuon trigger requirements. To avoid regions of 



rapidly changing trigger efficiency, we omit CMU (CMX) 
muon candidates with px < 2 (2.2) GeV/c. The random 
combinatoric backgrounds are suppressed by requiring 
the vector sum of the muon transverse momenta to be 
IPt^I > 4 GcV/c. The remaining pairs of muon tracks 
are then refit under the constraint that they come from 
the same three-dimensional (3D) space point, and are 
required to satisfy vertex fit quality criteria. The 3D 
decay length is given by — L ■ where 
L is the displacement vector from the primary to the 
dimuon vertex. The primary vertex is determined using 
a constrained vertex fit of all tracks in the event, ex- 
cluding the pair and other secondary decay tracks. 
For each B-candidate we estimate the proper decay time 
t = M mm L3d/|p mm |, where M MM is the invariant mass 
and p MM is the momentum vector of the dimuon sys- 
tem. Additional background is reduced by demanding 
L3D < 1.0 cm (to remove poorly reconstructed tracks), 
the uncertainty on L3D to be less than 150 /im, and 
2<j\ < A < 0.3 cm, where X — c ■ t and g\ is the to- 
tal uncertainty on A. There are 22459 (14305) dimuon 
candidates that fulfill all the above trigger and offline 
reconstruction requirements in the U-U (U-X) channel. 
At this stage, the data sample is dominated by random 
combinatoric background. 

We model the signal B° s d — ► decays using the 

pythia Monte Carlo (MC) program Q. The pythia 
events are passed through a full detector simulation and 
satisfy the same requirements as data. The px spectrum 
of the B-hadron in the MC is weighted to match the mea- 
surement from Ref. [lfj . The signal MC samples are used 
primarily for analysis optimization and for estimating the 
efficiency of selection requirements. 

For the final event selection we use the following four 
discriminating variables: M MM , A, the 3D opening angle 
AO between vectors y ^ and L, and the 5-candidate 
track isolation (I) 15]. We exploit long lifetime of B 
mesons to reject prompt combinatoric background with 
the decay length (A) requirement. Additional combina- 
toric background and partially reconstructed B hadrons 
are removed with the pointing angle requirement. Since 
&-quark fragmentation is hard, B hadrons carry most of 
the transverse momentum of the b quark, and thus are 
isolated. We use the variable I to enhance the heavy fla- 
vor content of the sample and also to reject partially re- 
constructed B hadrons, which are less isolated. Figure 
compares the distributions of these variables for data 
(which is background-dominated) to MC signal events. 
Based on the observed distributions, we apply two ad- 
ditional loose requirements, I > 0.5 and AO < 0.7 rad, 
to further reject background while maintaining 92% of 
signal efficiency. We collectively refer to all the selection 
requirements applied up to this point as the "baseline" 
requirements. In the data, 6242 (4908) events survive the 
baseline requirements in the U-U (U-X) sample. 
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FIG. 1: Distributions of the discriminating variables for data 
(dominated by background) shown in solid histogram and sig- 
nal — > fj, + events shown in dashed histogram. 
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FIG. 2: The likelihood ratio distribution for signal (dashed) 
from pythia MC and background (solid) from data sidebands 
in the U-U channel. Similar shapes are observed in the U-X 
channel. 



To enhance signal and background separation we con- 
struct a multivariate likelihood ratio based on the input 
variables: /, A9, and A probability P(A) = e ~ X,CTB ^ d ) , 
where T Bs{d) is the world average B s ^ lifetime. We use 
the -P(A) variable instead of A in constructing the likeli- 
hood ratio because the P(A) distribution is nearly flat, 
and therefore better behaved in the likelihood. We define 
the likelihood ratio to be 



L R = 



(i) 



where x\ = I, x% = A9, X3 = -P(A), and P s (6)(:Ci) is 
the probability that a signal (background) event has an 
observed X4. The probability distributions for the signal 
events are obtained from the signal MC and the back- 
ground distributions are taken from the data sidebands 
(defined below). The resulting Lr distributions for the 
signal and background events are shown in Fig. [3 

Although a subset of the data was used previously || , 
we adopt an optimization strategy that only uses events 
in the data sidebands in order to avoid biases in our 
choice of final selection criteria. The search window is 
defined by 5.169 < < 5.469 GeV/c 2 , corresponding 
to approximately ±6 times the two-track invariant mass 
resolution, which is estimated to be um « 24 MeV/c 2 . 
The sideband regions 4.669 < Af w < 5.169 GeV/c 2 and 
5.469 < M w < 5.969 GeV/c 2 are used to estimate the 
background in the signal region and for extensive analysis 
cross-checks. 

The Bg — * branching fraction is obtained by nor- 

malizing to the number of B + — » J/ipK + — > /i + /j,~K + 



decays collected by the same trigger. The B + — > J/ipK + 
mode is reconstructed us ing the same baseline require- 
ments as the signal mode 16] but including an additional 
Pt > 1 GeV/c requirement on the kaon candidate. The 
upper limit on the branching fraction is given by 



B(B° S 



/i+M - ) 90%c.i. 



N 9Q n % „ I 

iv sj a B + e B+ 1 



a B o 
k . B (B- 

Is 



-base 
t B 



■ B" 



J/^K+),(2) 



decays at 



where N^ % is the number of B® 
the 90% C.L. for N observed and Nb expected back- 
ground events. The value of N B °o /o is estimated using 
the Bayesian approach of Ref. [I3 assuming a flat prior 
for the B(B® — > /i + /i _ ) and incorporating Gaussian un- 
certainties into the limit. The number of reconstructed 
B + — > J/tpK + candidates, N B +, determined from the 
data using sideband subtraction, is 1785 ± 60 for the U- 
U and 696 ± 39 for the U-X channel. The parameter 
a B o (a B +) is the trigger acceptance and e b B o e (e b B + e ) is 
the efficiency of the baseline requirements for the sig- 
nal (normalization) mode. We apply the likelihood ratio 
requirement only to the signal mode and therefore the 
efficiency of the likelihood ratio e^J appears only in the 
denominator of Eq. |21 We use the fragmentation ratio 
f u /f s = 3.71±0.41 18]. The branching fraction B(B+ -> 
J/1PK+ -»■ uUjrK + ) = (5.88 ±0.26) x 10~ 5 is obtained 
from Ref. [17j. The expression for B(B ( j — > /i + /i _ ) is 
derived by replacing B® with B® and the fragmentation 
ratio with f u /fd = 1. 
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The expected number of background events is esti- 
mated by extrapolating the number of sideband events 
passing the baseline requirements to the signal window 
and then scaling by the expected rejection factor k for a 
given Lji requirement. The parameter k is determined 
from the background Lr distribution, which we generate 
by randomly sampling the A, A0, and / distributions 
from the data sidebands to improve statistical precision 
on k. The relative uncertainty on Nf, is ±15% (±19%) 
for the U-U (U-X) channel. The dominant contribution 
comes from the limited statistics of the input distribu- 
tions used to generate the background Lr distribution. 

We have cross-checked our background estimate pro- 
cedure using control samples from the data: like sign 
/i ± /i ± events, /i + \i~ events with A < 0, and a fake-muon 
enhanced sample in which we demand at least one 

muon candidate fail the muon quality requirements. We 
compare our background predictions to the number of 
events observed in the search window for a wide range of 
Lr requirements. No statistically significant discrepan- 
cies are observed. For example, requiring Lr >0.99 and 
summing over all control samples, we predict 3.8 (4.0) 
background events and observe 4 (3) in the U-U (U-X) 
channel. In addition we have considered background con- 
tributions from B® d — > h + h~ , where — n^ 1 or , in 
the signal window and determined that the contribution 
from those decays is negligible. 

The acceptance ratio ag+ / a B o obtained from the sig- 
nal MC is 0.297 ± 0.020 (0.191 ± 0.013) for the U-U (U- 
X) channel. The uncertainty includes contributions from 
systematic variations of the modeling of the B-hadron 
Pt spectrum and longitudinal beam profile, and from the 
statistics of the MC sample. 

The quantity e base includes the trigger and offline re- 
construction efficiencies. The trigger efficiency is deter- 
mined from inclusive data samples unbiased with respect 
to the triggers used here. The ratio of B + to B® trig- 
ger efficiencies is measured to be less than 0.1% away 
from unity. We evaluate the single track COT, SVX II, 
and muon efficiencies using a data sample of inclusive 
J /if) — » decays collected on single-muon triggers. 

The relevant double-track efficiencies are computed by 
convolution with B° -» ^+ ^ and B + -» J/ipK + MC 
events surviving the trigger requirements. The offline re- 
construction efficiency between signal and normalization 
mode also largely cancels in the ratio with the excep- 
tion of the kaon efficiency from the B + decay. Lastly, we 
obtain the efficiency of the remaining baseline require- 
ments from the signal MC and cross-check the results by 
comparing B + data and MC. Combining all effects, we 
find e b ^ e /e b ^ e = 0-920 ± 0.034 (0.915 ± 0.034) for the 

U-U (U-X) channel. The uncertainty is dominated by 
systematic uncertainties accounting for kinematic differ- 
ences between J /ip — » and B® d — > /i + /i~ decays. 

The efficiency of the likelihood requirement e 

d(? is es- 



timated from the signal MC. The efficiency varies from 
about 70% for Lr > 0.90 to 35% for L R > 0.99. 
We assign a relative systematic uncertainty of ±5% to 
both U-U and U-X channels based on comparisons of 
B + -> J/ipK+ MC and data samples. 

We optimize the analysis based on the a priori ex- 
pected 90% C.L. upper limit on B(B° sd -> fi + fJ,~). The 
expected limit for a given set of optimization require- 
ments is computed by summing the 90% C.L. limits over 
all possible observations N, weighted by the correspond- 
ing Poisson probability when expecting Nb. We scan over 
a range of Lr requirements and determine the optimal 
value to be Lr > 0.99. With the optimized selection re- 
quirements, the expected number of background events is 
the same in the B° s and B d search windows. Within the 
signal region of ±60 M eV/ c 2 (±2.5<7m) about the world 
average B° or B° d mass H3, N b is 0.81 ±0.12 [0.66±0.13], 
and the B® single-event-sensitivity is (1.0 ± 0.2) x 10~ 7 
[(1.5 ± 0.3) x 10~ 7 ] for the U-U [U-X] channel. Using 
these criteria we observe no events inside either the B® 
or B d signal box as shown in Fig. [3J The one candidate 
event that survived in the previous analysis has an 
Lr value of « 0.8, and fail the selection requirements 
here. Using Eq. [21 and combining the U-U and U-X 
channels taking into account the correlated uncertainties, 
we derive 90% (95%) C.L. limits of B(B° -> < 
1.5 x 10~ 7 (2.0 x 10~ 7 ) and B{B G d -> ^ + M") < 3.9 x 10~ 8 
(5.1 x 10 -8 ). The new limits improve the previous lim- 
its 0, 3 by a factor of two and can be used to reduce the 
allowed parameter space of a broad spectrum of SUSY 
models |Ellli 
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FIG. 3: The fi + n~ invariant mass distribution versus likeli- 
hood ratio for events satisfying baseline requirements for U-U 
(solid triangle) and U-X (open circle) channels. The B° (solid 
box) and (dashed box) signal regions are also shown. The 
one candidate event observed in the previous analysis is high- 
lighted with a star symbol (at Lr ai0.8). 
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